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RADIATIVE INTERACTIONS BETWEEN ABSORBING-EMITTING AND 
FLOWING MEDIA WITH INTERNAL ENERGY GENERATION 
by John R. Howell 
Lewis Research Center 

SUMMARY 

An analytical study of the local energy flux and emissive power in two adjacent re- 
gions is given. Each region is assigned a radiation absorption coefficient, mass flow 
rate, heat capacity, and volumetric energy source. The effect of these parameters and 
of the relative thickness of the regions is examined. The regions are assumed to be in- 
finite, parallel layers of finite thickness in one case, and coaxial cylinders of finite radius 
in another. The problem of determining the local one-dimensional emissive power and 
radiant energy flux is solved by using a second-order diffusion solution with emissive 
power jump boundary conditions where applicable, and in the parameter ranges where 
such an approach is invalid a Monte Carlo transport analysis is used. The diffusion solu- 
tion for the cases examined is accurate for optical thickness greater than about 1. The 
regions are assumed to be gray and have constant properties. A sample problem which 
demonstrates the use of the derived equations is presented. 


INTRODUCTION 

Radiative interchanges between adjacent absorbing- emitting media with different 
properties and internal energy sources are being examined because of applications to a 
number of practical systems. These systems include radiative shielding by injection of 
seeding material, radiative transfer in the gaseous-core nuclear propulsion system, and 
hypersonic shock phenomena. For many such problems, radiation is the major mecha- 
nism for energy exchange. The present work provides a simplified solution to radiative 
exchange problems where two distinct media are present. From this the engineer may 
estimate the importance of individual parameters on radiative transfer for his specific 
problem. 

A diffusion solution including second- order terms and the effect of temperature slips 



at region boundaries (ref. 1) is extended for use in two geometries. This type of ap- 
proach leads to closed solutions in the form of algebraic equations and provides good ac- 
curacy over a very wide range of parameters (a much wider range than is generally con- 
ceded). Inclusion of the second-order and slip effects extends the results to most regions 
common to the problems mentioned previously. Certain regions do remain, however, 
where diffusion solutions are of doubtful validity or fail altogether. For this reason, a 
transport solution using a Monte Carlo technique (refs. 2 to 4) is utilized for determining 
where the diffusion solutions deviate appreciably from the more exact results. 

All solutions, for both generality and simplicity, have the implicit assumptions of 
slug flow profiles, local thermodynamic equilibrium, no scattering, constant properties, 
and gray gases. In addition, the change in temperature with distance parallel to the flow 
direction is considered constant. Any or all of these may be very poor assumptions for 
a given problem, yet to remove any one of them gives a resulting loss of applicability to 
some problems otherwise covered. The assumption of a gray gas is a case in point. 
Various systems being considered, notably those involved in gaseous-core nuclear pro- 
pulsion, might use small particle seeds to attain radiative absorption and can be analyzed 
on a gray gas basis, with some justification, so long as the wavelength of the radiation is 
not large compared with the particle size. Just as large a variety of problems, however, 
involve real gases where a gray gas assumption is simply incorrect. To allow solution of 
any of this latter class would probably mean a specialized approach to the exclusion of 
many related problems. For this reason, the assumptions are included, with the caution 
that use of the equations presented herein implies a responsible look by the user at the 
assumptions involved. 


SYMBOLS 


A area 

C heat capacity 

Jr 

D ratio of plane region thicknesses, Lj/Lg 
D c ratio of cylinder diameters, Dj/Dg 
E radiant energy per unit time 

4 

e emissive power, aT 

e^ spectral blackbody energy distribution (spectral emissive power) 

I mean square absorption coefficient defined by eq. (A12) 

L dimension in z direction 

£ region thickness 
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Q volumetric energy generation rate 

q energy flux 

R dimensionless radius, r/Rg 

r radial cylindrical coordinate 

T absolute temperature 

u velocity of flow 

V volume 

x Cartesian coordinate in flow direction 

y Cartesian coordinate perpendicular to x and z 

Z dimensionless coordinate, z/Lg 

z Cartesian coordinate normal to bounding plane 

r total energy sources per unit volume, -pC u(dT/dx) + Q 

e emissivity of bounding surfaces 

k gas absorption coefficient 

Rosseland mean gas absorption coefficient, eq. (A10) 

k mean absorption coefficient, eq. (All) 

s 

p density 

a Stefan- Boltzmann constant 

r optical thickness, k!L 

<p(P) <p(Z) evaluated at Z = D 

cp( R) dimensionless local emissive power for concentric cylinders, [e(R) - e w ]/r^Rj 

<p( Z) dimensionless local emissive power for parallel planes, [e(Z) - e^/r^L^ 

xp(D) i//(Z) evaluated at Z = D 

t//(R) dimensionless energy flux in concentric cylinders, -q/r^Rj 
i//(Z) dimensionless energy flux in parallel planes, -q/r^Lj 
S2 function defined by eq. (A4) 

a) solid angle 

Subscripts: 

D evaluated at Z = D 
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g gas 

g-g dimensionless gas interface slip condition 
g-w dimensionless wall slip condition 
g, w evaluated in gas very close to wall 
w bounding wall black emissive power 
x in x- direction 

z in z -direction 

X spectral value 

0 origin for Taylor series expansion 

1 region 1 

2 region 2 

ANALYSIS 

The diffusion equations for the radiative 
interchange between two adjacent absorbing- 
emitting regions are derived in the appendix, 
and the resulting equations for the local 
emissive power and radiant flux are summa- 
rized in this section. Two geometries are 
examined. 

The first arrangement consists of two 
parallel planes of finite thickness as shown 
in figure 1 (a) . The interface between the re- 
gions is assumed to be an imaginary plane 
across which no mixing occurs and which 
offers no resistance to radiative energy 
transfer. The outer boundary of region 2 
(z = Lg) is assumed to be a solid gray wall 
of emissivity e . The inner boundary of re- 
gion 1 (z = 0) is taken to be a plane across 
which there is no net radiative flux. The two 
planes examined here may then be viewed as 
one-half of a symmetrical system made up 
of a central or core region flanked by two 
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identical regions with properties differing from the core. 

The second geometry is shown in figure 1(b) and consists of two concentric regions 
of absorbing- emitting media bounded on the outer surface (r = rg) by a solid boundary of 
emissivity e . The interface between the regions is an imaginary cylindrical surface 
assumed to have no radial mixing occurring across it and offering no resistance to radi- 
ative transfer. 


Diffusion Solution 

Parallel plane geometry . - In the appendix, the second-order diffusion solution as 
proposed by Deissler (ref. 1) is used to determine the one- dimensional emissive power 
distribution for the two cases. First, the parallel layer geometry (fig. 1(a)) with slug 
flow traveling parallel to the boundaries and with constant and uniform internal energy 
generation in each region is examined. Individual absorption coefficients for radiation, 
k j and Kg, are present in regions 1 and 2, respectively. If the dimensionless parame- 
ters are defined as 


e_(Z) - e_ 
<p(Z)=-£ £ 

r l L l 



t 2 = *2^2 " L l^ 

T = -oC u — + Q 
dx 

then the emissive powers in regions 1 and 2 at point Z are given by equations (A29) and 
(A26) of the appendix, respectively, which are 
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0 < Z < D (A29) 


3ti p 

<Pl ( z ) “ 6 <P g _ w = ? 2 (°) - 6<p g-w + — 4 t° 

8D 2 


- z2 > + e * e -g 


and 


3r 2 (r 2 /r 1 )(i - Z 2 ) 3t 2 

<P 2 (Z) ~ 6< ?g-w = ~ + — I 1 " — I d - z ) 

2 g w 8D(1 - D) 4(1 - D) V I\ 


D< Z < 1 


(A26) 


where the slip values across the gas-wall and gas-gas interfaces (5<p g _ w and 6< P g _ g> 
respectively) are given by 


and 


6<p g-w 


l_l\/l-D r 2 \ , 3 ( r 2/ r iH 1 - D > 

,e 2/^ D T 1 y 8 Dt 2 


(A27) 


6<p 


g-g 


_3_ 

8D 



1 - D 

t 2 



(A28a) 


Equation (A29) shows that the local dimensionless emissive power of the gas in re- 
gion 1 depends on the emissive power in region 2 evaluated at the boundary between the 
regions, the relative thickness D, the optical thickness Tj, and 6<Pg_g* The last term 
accounts for the discontinuity in emissive power that occurs at the boundary between 
absorbing- emitting regions with internal sources (refs. 5 and 6). 

The local emissive power in region 2 depends on the optical thickness r 2 in region 2, 
the relative thickness D of the layers, and the ratio of sources in the two regions r 2 /Tj. 
In addition, there is a dependence on the emissive power discontinuity at the bounding 
surface 6<P g _ w , which, in turn, depends on the same parameters and the wall emissiv- 
ity e. 

The local radiative flux in a direction normal to the slab boundaries is given by equa- 
tions (A36) and (A3 5) in the appendix as 


(Z) = Z 0 < Z < D 
1 D 


(A36) 
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D < Z < 1 


^ 2 (z) = 1 + — ” 

r i 



(A3 5) 


where <//(Z) is a nondimensional radiative flux given by 

i^(z) = -S@L 

r i L i 

Concentri c cylindrical geometry . - It is also shown by equations (A40) to (A43) in the 
appendix that, for the geometric model shown in figure 1(b), the emissive power distribu- 
tion radially is given by 


3 t ^2 ^ j ^2 

(0j<r) - t<p g . w - <p 2 ( d c > - + 5?, g-g 

®c 


0 < R < D c (A40) 


and 


<p 2 (R) - 6 <p g _ w 


3t 2 r 2 


1 - R 


16 T t D c (l - D c ) 


3TnD / I^oX 

(1 - — “ ]ln R 

8(1 - D c ) y TjJ 


D,SRSl (A41) 


where 


%-w 



1\ A 2 1 - D c 

2 A r l D c 



3 1 - D c / r 2 1 •> D c 

16 r 2 \Tj D c 


and 


sr-er 
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g-g 32D, 



(A42) 


(A43) 


The dimensionless variables have the same physical interpretation as before, but now 
they have slightly different forms in some cases to reflect the difference in geometry. 
The variables are 
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e (R) - e 

<P( R)=- g 

r i r l 


D c=- 

r 2 


T 1 = K l r l 


t 2 = K 2^ r 2 ‘ r l ) 


R = 


The local radial radiative flux is given by equations (A44) and (A45) as 


1/^01)=-^- 0<R<D 

2D„ c 


(A44) 


and 


d r 2 R 2 - D 2 

^ 2 (R) = — + — ^ 

* 2R Tj 2RD 


~ R “ 1 


(A45) 


where 


i//(R) = --SM 

r iri 


Monte Carlo Transport Solution 

For comparison with the diffusion solutions presented in the preceding sections, a 
transport analysis was made by means of Monte Carlo techniques. For the parallel plane 
geometry, the essential equations and concepts are available in reference 2. For the 
cylindrical case, similar information is available in reference 3, and a general discus- 
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sion of Monte Carlo approaches to radiative transfer problems is presented in refer- 
ence 4. 

Because of the need to define volume elements for determining local values of gas 
emissive power, it is not feasible to determine the gas emissive power exactly at the 
gas -wall boundary. Therefore, it is assumed in the Monte Carlo solutions presented 
herein that the emissive power slip at the wall could be adequately represented by extrap- 
olation of the radial gas emissive power curve to the wall from the three calculated 
<p 2 (Z) points in the gas adjacent to the wall. 

Because the Monte Carlo technique involves the sampling of a model of the radiant 
energy exchange process by statistical means, the answers obtained contain statistical 
fluctuations. Some estimate of the probable error in the results is desirable. The 95- 
percent confidence limits on the Monte Carlo points are therefore used to indicate that 
the value of the point has a probability of 95 percent of lying within the range shown. The 
limits did not usually exceed the size of the symbols used to plot the points. Where they 
did, the limits are indicated by vertical bars. 



Figure Z - Dimensionless gas emissive power distribution in 
gray planes. Ratio of energy sources, 1; ratio of thicknesses, 
0.5; optical thickness in region 1 same as in region Z 


RESULTS 

Comparison is now made of the dif- 
fusion and transport solutions to deter- 
mine the range of parameters for which 
the diffusion solution is valid. It is ex- 
pected from the assumptions inherent in 
the diffusion approach that diffusion re- 
sults would be most in error for small 
values of r and in regions within one 
mean free path of a boundary. 

Parallel Planes 

In figure 2, both regions are given 
identical properties and equal widths 
(D = 0. 5). The effect of optical thickness 
on the emissive power distribution is then 
demonstrated. Accuracy of the diffusion 
solution on a percentage difference basis 
is seen to be excellent above r ~ 2 and 




Figure 3. - Effect of optical thickness on radiation slip. 
Ratio of thicknesses, 0.5; optical thickness in region 1 
same as in region 2. 



Ratio of energy sources, r 2 /I" j 


Figure 4. - Effect of energy source ratio on emissive power slip at wall 
in parallel plane regions. Optical thickness in region 1, 2; optical 
thickness in region 2, 2; thickness ratio, 0.5. 


acceptable at even smaller values of t. 

In figure 3, the emissive power 
slip at the solid boundary is shown as 
a function of gas optical thickness in 
region 2, with the ratio of energy 
source strengths as a parameter. 

Again, accuracy is quite good for the 
diffusion solution at t^ = 2 and ap- 
pears to diverge from the Monte Carlo 
results below this value. Even at 
Tg = 1, however, the errors between 
the results are small. It should be 
noted that the equation for 6<pg_ w by 
the diffusion approach does not depend 
on the energy reaching the bounding 
surface from region 1. Yet, for a 
very small Tg, it would be expected 
that 6<pg_ w would depend chiefly on 
that energy, so that certain combina- 
tions of parameters might lead to more 

error in 6 cp than is implied in fig- 

g w 

ure 3. 

In figures 3, 5, 8, and 10, the 
source strength ratio I^/T ^ appears 
as a parameter. The definition of r 
is the net volumetric rate of energy 
accumulation due to both internal en- 
ergy generation and change in the en- 
thalpy of the media. If the flow con- 
ditions of the media are such that the 
enthalpy term is greater than the en- 
ergy generation term, then it is quite 
possible for T to have negative values 
in one or both regions. This simply 
means a net volumetric energy sink is 
present in the region. Both positive 
and negative ratios of source strengths 
are therefore shown. 

The effect of the energy source 
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Figure 5. - Effect of energy source ratio on emissive power. 
Optical thickness in region 1, 2; optical thickness in 
region 2, 2; thickness ratio, 0. 5. 


ratio of the regions Tj/rg on the emissive 
power slip at the wall is demonstrated in 
figure 4. Three results are noted. First, 
the difference between the diffusion and Monte 
Carlo results increases as the absolute mag- 
nitude of ^ 2 /r 1 increases. Second, the 
absolute magnitude of the slip itself increases. 
Third, the percentage difference between the 
two solutions does not change markedly as 
the source ratio becomes large. 

Figure 5 shows the emissive power dis- 
tributions with the slip value subtracted out 
to avoid a greatly expanded scale. This is 
valid because both methods were shown to 
predict accurately the wall slip in emissive 
power. Therefore, it is evident that large 
energy source ratios give larger differences 
between the solutions. However, if the slip 
values are added back in, the percentage 
difference between the solutions is not a func- 
tion of source ratio, as far as can be deter- 
mined. This trend continues out to source 
ratios of plus and minus 10^ (for which re- 
sults are not presented). This subtraction 
of the slip values also accounts for the ap- 
parent dependence of the percentage differ- 
ence on the sign of which does not 

appear if the slip is added to the results. 

It is interesting to note the emissive 
power jumps at the region boundary. The 
magnitude of these jumps increases with in- 
creasing magnitude of source strength ratios. 
This is in accord with the results of refer- 
ence 6. 

The remaining parameter to be examined 
is the slab thickness ratio D. The effect on 
c p (at two positions) of varying this parameter 
while holding all other conditions in the re- 
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gions constant is shown in figure 6 for one 
case. The absorption coefficients and 
*2 are constant, so that Tj and Tg are now 
functions of D. It is seen that for intermedi- 
ate D values good agreement between solu- 
tions is reached. This is expected because 
the diffusion solution, even including the emis- 
sive power jump at the bounding surface, 
becomes inaccurate within a radiation mean 
free path of a boundary. The best agreement 
should occur for Tg = which is the case at 
D = 1/6. It is difficult from the results given 
to determine whether this is the case. It is 
interesting that for values of D ~ 1, the dif- 
fusion solution will give a slip prediction 
(eq. (A27)) based only on e and Tg, while in 
actuality the slip for D ~ 1 should be a func- 
tion also of conditions in region 1 as the opti- 
cal thickness of region 2 becomes small. 

This is why <p<£L — 1) for the diffusion solu- 
tion begins to diverge from the Monte Carlo 
analysis as D approaches 1 in figure 6. 

For the parallel plane geometry, then, 
the greatest inaccuracies in the diffusion re- 
sults occur for small r, or for D near 0 or 1. This would be predicted a priori from 
knowledge of the assumptions inherent in the diffusion approach. What is surprising is 
that the errors are of the order of a few percent even for r « 1 and for extremal values 
of D. 

In figures 7 to 11, an analogous set of results is presented for the case of concentric 
regions as shown in figure 1(b). The conclusions remain the same as for parallel regions, 
in that the accuracy of the diffusion solution is very good for all the cases examined. Ac- 
curacy for like values of parameters is slightly inferior to the results for parallel regions. 


Thickness ratio, D 

Figure 6. - Effect of plane region thickness on emissive power. 
Energy source ratio, 1; optical thickness in region 1, 10D, 
where D is thickness ratio; optical thickness in region 2, 
2(1 - D). 


Sample Problem 

To illustrate the use of the diffusion solution in a more practical calculation, the 
radial temperature distribution in a coaxial- flow gaseous-core nuclear propulsion system 
will be determined. Such a device consists of a core of fissionable material in the gase- 
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Dimensionless emissive power, 9-69 



Position in medium, R 

Figure 7 . - Dimensionless gas emissive power distribution 
in gray coaxial regions. Ratio of energy sources, 1 ; 
optical thickness in region 1 same as in region 2. 


Diffusion Monte Carlo Energy source 
results results ratio, 

Vl 



2 1 1 L I I I 

0 2 4 6 8 10 

Optical thickness in region 2 , T2 

Figure 8. - Effect of optical thickness on radiation slip at 
wall bounding coaxial regions. Ratio of cylinder diam- 
eters, 0 . 5 ; optical thickness in region 1 same as in 
region 2. 


i 



-4 0 4 8 12 

Ratio of energy sources, 

Figure 9 , - Effect of energy source ratio on emissive power slip at wail 
in concentric flow regions. Optical thickness in region 1 , 2 ; optical 
thickness in region 2, 2; ratio of cylinder diameters, 0. 5. 
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Position in medium, R 

Figure 10. - Effect of energy source ratio on emissive power 
in coaxial regions. Optical thickness in region 1, 2; 
optical thickness in region 2, 2; ratio of cylinder diam- 
eters, 0.5. 
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Diameter ratio of cylinders, D c 

Figure 11. - Effect of diameter ratio on emissive 
power. Energy source ratio, 1; optical thick- 
ness in region 1, 10D C , where D c is ratio of 

cylinder diameters; optical thickness in 
region 2, 2(1 - D c ). 


ous state, generating heat internally by the fission process. The heat is transferred to 
an annulus of flowing propellant, which is then expanded through a nozzle to provide 
thrust. 

The following are the conditions for study based upon some more or less typical con- 
ditions (ref. 7). The reactor core has a diameter of 3. 2 feet, and the outer diameter of 
the propellant annulus is 8 feet. Then 

D = 3 ' 2 / 2 =0.4 

c 8/2 

4 

A fission material of optical thickness = 10 and a propellant seeded to Tg = 7 are 
used. Also, the propellant flow is adjusted so that 90 percent of the radiant heat gener- 
ated in the core is absorbed in the propellant. Then 

0. 9 » -r 2 ir (l - D^) 


or 



-0. 172 
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200xl0 3 



figure 12. - Radial temperature dis- 
tribution in sample gas-core reactor. 

Finally, an emissivity € = 0. 8 can be assigned to the bounding surface of the system as 
a reasonable choice for most metals at high temperature, and the emissive power distri- 
bution in the core and propellant regions can be computed by use of equations (A40) to 
(A43). To find the temperature distribution, the volumetric heat generation rate Tj in 
region 1 must be specified (eq. (A44)). In the real case, a number such as = 1. 4xl0 5 
Btu per cubic foot-second is in the region of interest. With this value, the computed 
temperature distribution is given in figure 12. 

In any calculation used in practical design, account must be taken of the property 
variations with temperature. Allowing k to vary with temperature in more sophisticated 
calculations brings the predicted temperatures in the core down to a maximum of about 
100 000° R (ref. 7), which is a considerable change from the prediction of 172 000° R in 
this investigation. 

The calculation performed herein can be done without a computer and gives an indica- 
tion of the problems to be expected, that is, the very steep temperature gradient in the 
fuel near the region boundary and the very high temperatures at the center of the fueled 
region. The results were calculated by using equations (A40) to (A43) without the neces- 
sity of resorting to more complex transport analysis. 


CONCLUSIONS 

The second-order diffusion solution incorporating jump conditions at region bound- 
aries as proposed by Deissler was extended to problems involving two distinct media, 
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and its region of validity was determined by comparison with a Monte Carlo transport 
solution. 

The emissive power jump that occurs between two absorbing- emitting regions can 
become important over certain ranges of parameters, and should be examined in any situ- 
ations where it may exist. The diffusion solution yields accurate answers when applied to 
problems involving two regions undergoing radiative interchanges of energy, and it can be 
used over a wider range of parameters than might have been expected. For the problems 
examined, the diffusion results were within a few percent of the more accurate transport 
solution so long as the optical thickness in each region was greater than about 1. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 8, 1966, 

129-01-09-09-22. 



APPENDIX - DIFFUSION SOLUTION 


As shown by Deissler (ref. 1), the emissive power of a gas near a point x Q , y Q , z Q 
can be expanded in a Taylor series to give 



By writing equation (A2) in spherical coordinates, applying the equation for emission of 
energy from an element which arrives at the element of volume around Xq, y^, as 


dE. = 4 k e dV — e ^ 
A A A 4jf 


and using these to determine the net flux per unit area per unit frequency interval pass- 
ing through element dA in direction z, Deissler shows that the net flux is given by 



where 


n(n, v, s) 


♦ <-l) T - s ] [l * <-l) s ] 

n! r^ n “ v + 

r ( 

V - s + l'l 

2 J 

r < 

(■,*■) 

(n - v)! (v - s)! s!T| 

ft 

+ 4' 
2 / 



(A4) 


18 



The r in this equation is the gamma function and is not to be confused with the source 
strength used elsewhere in this report. For a wall bounding the gas at some positive 
value of z, the energy jump at the wall is given by 



To this point it has only been assumed that the absorption coefficient of the gas does not 
vary over a mean free path and that the flux solution applies throughout the gas, although 
it is strictly correct only one mean free path, or more, away from a bounding surface. 

If terms with derivatives of higher order than 2 are neglected, equation (A3) reduces 
to 




_4_^X 
3k^ 3z 


(A6) 


and equation (A 5) becomes 



Writing Planck' s spectral distribution function in terms of black emissive power 
rather than in terms of temperature and integrating equations (A6) and (A7) over 
wavelength interval AX gives 


4 

e_ = ctT* 
g 
the 


3e 


g 


3 k r , ax 9z 


(A8) 


and 


19 




To this point, the development is essentially that presented by Deissler (ref. 1) as 
extended in the discussion appended to reference 1; it differs from the usual diffusion 
solution chiefly in the retention of second-order terms and in the inclusion of the emissive 
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power jump (eq. (A9)) at the solid boundary. However, the changes in the results due to 
these two factors are significant, as is shown by the comparisons in references 1 to 3. 
These references are for gray gases. When wavelength effects become important, and 
when is a function of temperature, agreement with exact solutions may become less 
satisfactory, as shown in reference 4. 

For the problem at hand, derivation of one more condition is necessary. At the 
boundary between two regions with different absorption coefficients (fig. 1), the total 
radiant energy per unit frequency interval passing from region 1 to region 2 through an 
area dA on the boundary of the regions is 



and from region 2 to region 1 it is 



The net radiant flux per unit frequency interval in the z-direction is then 



Neglecting terms of order higher than 2 and integrating over all wavelengths as before 
give 
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r 



For a gray gas, equations (A10) to (A12) give 


-i- = — and 1 = 0 
k R k s 


Plane Layers 


Distribution of emis sive power . - If attention is restricted to the model in figure 1(a) 
and if a gray gas and constant regional properties are assumed, the terms involving de- 
rivatives in x and y are zero, and equations (A8), (A9), and (A14) give, respectively, 



(A15) 


(A16) 
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(e _ e ) =-q + J- _L /_!] +J_ /_i 

g,2 g.lg.g z >° I 3 ^^as/g Kj \3z 


1 1 

+ — 


2 * 2 1 V z 7l *2 V z 72 


(A17) 


In both layers 1 and 2, the net heat flow out of a unit volume must be equal to the volumet- 
ric heat generation rate Q minus the rate of change in enthalpy due to flow in the volume. 
If steady conditions are assumed and flow occurs only in the x-direction, then, for layer 2, 


• < l Z ,2= Q 2-(' 5C p u x-^) L 1 szsL 2 




3k 2 0 z 2 


(A 18) 


Putting equation (A18) into dimensionless form gives 


3 e ^ r 2^2 ^2 

g _ z z z D < Z < 1 

a7 2 4(1 - D) 


(A19) 


where 


Z T 2-'‘2< L 2- L l>; 


In region 1, a similar analysis gives 


a z e„ 3T.r-.Ln 
g= 1 1 2 0<Z<D 
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Integration of equation (A20) with respect to Z gives 


Be 3 t . r ^ LnZ 

— £ = — LJ_A_ + Cl 

3Z 4D 1 


where is a constant. However, at Z = 0, (Be^/BZ) = 0 because of the zero flux con- 
dition for the model chosen, so that C ^ = 0 and 



az 


3r i r l L 2 Z 

4D 


0< Z < D 


(A21) 


Now, integrating equation (A19) with respect to Z gives 


9e g _ 3T 2 r 2 L 2 Z 
3Z 4(1 - D) 


+ Cr 


D < Z < 1 


(A22) 


To evaluate the constant C 2 , it is apparent from equation (A15) that if the radiant flux 
varies continuously between two regions with different absorption coefficients, then the 
gradient in emissive power at the boundary between the two regions is related by 


% 




or, in dimensionless form, 


D_ / 9e g\ = 1 - D / 9e g \ 

T 1 \ aZ /l,D t 2 \ 9Z /2,D 


(A23a) 


(A23b) 


Then, combining equation (A22) evaluated at Z = D with equations (A23) and (A21) evalu- 
ated at Z = D gives 


D 

C 2 L 2 ( r i ~ r 2> 

& 4 1 - D “ A A 


and 
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3e. 


D < Z < 1 


(A24) 


3t 2 L 2 I V Z " 3t 2 D 


3Z 


4(1 - D) 


4 1 - D 


L 2 r l 


Integrating once again gives 


e 


g 


3T 2 L 2 r 2 
4(1 - D) 



3t 2 D 
4 1 - D 


L 2 r l Z + C 3 


D< Z < 1 


(A25) 


The fact that at Z = 1 the value of eg is greater than the emissive power of the wall 

and will have a value of e_ denoted as e permits evaluation of the constant Co, 

S w ° 

and equation (A25) can be written as 


<?2< Z) ' 6<p g-w 


3 T 2(r 2 /r 1 )( 1 - z 2 ) 3 t 2 / r,\ 

8D(1 - D) 4(1 - D) y r J 


D< Z < 1 


(A26) 


where 


e„(Z) - e m 
<P( Z)=-8 

r l L l 


e - e 
6 cp 

^g-W -p T 

r i L l 


and e is the black emissive power of the bounding surface. The value of 6 <P~,„ is 

W &“ ” 

found directly from equation (A16) by noting that the net flux at the wall q z can be 
obtained from equations (A 15) and (A24). Substituting equations (A 15), (A24), and (A 19) 
into (A 16) and nondimensionalizing give 


6<p g-w 


'l l\ /l - D r 2 ( A t 3 ( r 2 / r i)( 1 ~ p ) 

S " 2/ \ D ) 8Dt 2 


(A27) 


The value of e. , n can be found in terms of e « n , which is the emissive power 
of the gas in region 2 evaluated at point D, and of the discontinuity or jump in emissive 
power between the two regions. The jump in emissive power at point Z = D can be 
found from equation (A17) after the emissive power distribution in region 2 is found. 
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Substituting both the value for q z D and also the appropriate derivatives into 

equation (A17) and nondimensionalizing give 


6<p 

g-g « d \ t 1 T v 


(A28a) 


or, in dimensional form 


g, l,D" e g,2,D 8 I 


3/ r i 


1 k 2 / 


(A28b) 


If the regions differ in the magnitudes of either their source strengths or their absorption 
coefficients, or both, a slip will occur in the gas emissive power at the boundary. The 
slip will only disappear in the special case Tj/kj = T^/k^ or if there is no internal 
generation in either region. 

Now the emissive power distribution in region 1 may be found by integrating equa- 
tion (A21) and using the slip condition of equation (A28). This gives 


3t, 

?1<Z) - 6 <Pg_ w = <P 2 ( D > - 6 V W + — “ (D* ~ Z‘) + 0 < Z < D 

8D 2 


(A29) 


The equations for <p in the plane slab geometry are then 


3t 


<pi(z) ' 5 V w = ?2< D > - 5 Vw + ~ (° 2 ~ z2 ) + 6 Vg 0 - z - D < A3 °) 

8D 2 


cp 2 (Z) - 6<p 


3 T 2 ( r 2 / r i) 

S- w " 8D(1 - D) 


(i-zV 


3 tc 


4(1 - D) 


(l - — ) (1 - Z) D < Z < 1 


g V 2/ \ r l D J 8D7 2 


(A31) 


(A3 2) 


6 cp 


3 /D 1 - D l 2' 


8D It, 


t 2 r l J 


(A33) 
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Energy fl ux distribution . - Combining equations (A23a) and (A26) gives 


or 


_ 4(1 - D) / 3e g\ 
3 t 2 \0Z/. 

_ 4(1 - D) 
w 3t 2 

_ 4(1 — D) r L 

3T 2 ( r 2 / r i) 

3t 2 

4D(1 - D) 


3 Tc 


Z + 


4(1 - D) 


1 


(A3 4) 


q 2 (z) r 2 ^ 


ipJz) = — = _e ± z + l — - 

" m T t~i tv -p 

A 1 


r l L l r l D 


= 1 + 


^(H dszs 


(A3 5) 


A similar combination of equation (A23a) evaluated in region 1 and equation (A21) gives 


^j(Z) = - 0 < Z < D 


(A3 6) 


Coaxial Flow Bounded by Cylinder 

Emissive po wer distribution . - For the system shown in figure 1(b), the energy equa- 
tion may be written as 

±Ai_4^i\ = _r (A3 7) 

3/c r dr y 9r J 


and (from eq. (A9)) the emissive power jump at the wall can be written as 


The emissive power jump at the boundary between regions is 



(A3 9) 


By following a derivation similar to that for plane layers, and noting that r in the equa- 

2 2 ^ 

tions must be replaced by (z + y ) in the evaluation of derivatives, it is evident that 


3 T j D 2 - R 2 


*’l< R )-%-w=*’2< I> c)- 6 »’g-w +5 ' , ’g-g + 7r-S- 0sEsD c < A40 > 


W D 2 


3t 2 r 2 1 - R 2 3t 2 D c 




c> S* 1 - D c> 



- — In R D < R < 1 (A41) 


where 


6cp 


g-w 


= 1 4 . i\ (hii^ +D V±i^ (h 1 ^ - 

2{e 2/ yTj D c C J 16 r, ^ D c C J 


6 ( P<r-<r = 


g_g 32D„ 


d ,i-d / r, 

Sl - 1 9 |4 — - - 1 ] 

T 1 3 t 2 


where 


(A42) 


(A43) 


eJZ) - e 


„ C ^ 1 X 7 1 

R = — ; cp{ Z) = -M ; D c = — ; t 1 = k^; r 2 = K 2 (r 2 - r^ 

r 2 r i r l r 2 
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Energy flux distribution . - By analysis similar to that for the plane geometry, the 
energy flux distribution is found to be 


^(R) = -5_ 0 < R < D„ (A44) 

2D c 


D r 2 R - 

i//«(R) = _£ + _£ S D < R < 1 (A45) 

* 2R T 1 2RD c 
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